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A  MODEL  OF  ACOUSTIC  BACKSCATTER 
FROM  ARCTIC  SEA  ICE 


1 .  INTRODUCTION 

1.1  Backscatter  levels.  Acoustic  backscatter  levels  from  Arctic 
sea-ice  have  been  observed1- **  to  be  much  stronger  than  observed  levels 
in  the  open  ocean.  Mellen5  demonstrated  that  observed  levels  were 
inconsistent  with  the  predictions  of  the  perturbation-type  spectral 
scattering  theories  from  a  rough  pressure  release  surface,  published 
by  Marsh”  and  later  by  Bass  and  Fuks7.  The  major  failing  of  the 
perturbation  models  is  to  neglect  the  effect  of  large-scale  surface 
features. 


1.2  Composite  surface  method.  Kuryanov8  and  Brown9  have  devised 
the  composite  surface  method  to  combine  the  effects  of  large-scale 
surface  features  and  small-scale  roughness.  In  Brown's  version  of  the 
composite  surface  theory,  the  surface  roughness  spectrum  is  divided, 
at  a  particular  wavenumber,  into  low-  and  high-wave number  portions. 
The  low-wavenuraber  portion  is  used  to  cal  alate  the  variance  of  a 
Gaussian  distribution  for  the  slope  of  the  surface.  The  high- 
wavenumber  portion  is  used  in  a  perturbation-type  model  to  calculate 
diffractive  scattering  strength  for  given  in-coming  and  out-going 
directions.  These  in-coming  and  out-going  directions  are  defined 
relative  to  the  scattering  surface,  and  hence  the  local  scattering 
strength  depends  on  the  local  slope  of  the  scattering  surface.  The 
mean  backscattering  strength  predicted  by  the  composite  surface  tech¬ 
nique  is  the  expected  value  of  the  local  scattering  strength  with 
respect  to  the  Gaussian  distribution  of  slope. 

1.3  Triangular  ridge  model.  Unfortunately,  the  slope  distribu¬ 
tion  of  Arctic  ice  is  not  Gaussian  in  nature.  The  composite  surface 
theory  with  a  Gaussian  distribution  of  slope,  while  giving  some 
improvement,  still  fails  to  predict  observed  backscatter  levels.  In 
the  next  section,  a  composite  surface  theory  based  on  scattering  from 
rough  triangular  ridges  will  be  derived.  The  triangular  ridge  model 
gives  good  agreement  with  backscatter  data  from  Arctic  sea  ice. 

2.  THE  SLOPE  MODEL 

2.1  Ridge  keel  structure.  The  underside  surface  of  Arctic  sea- 
ice  consists  of  relatively  smooth  stretches  of  ice  interrupted  by 
large  linear  ridge-keel  structures.  These  keels  are  piles  of  ice 
rubble,  with  steep  flanks.  In  the  triangular-ridge  composite  model 
proposed  here,  the  ice-ridge  keels  are  modeled  as  randomly  oriented 
prisms  of  triangular  cross  section,  whose  flanks  have  some  fixed 
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slope  angle.  Superimposed  on  Che  triangular  structure  is  small-scale 
roughness,  characterized  by  the  large-vavenumber  portion  of  the  sur¬ 
face  roughness  spectrum.  This  small-scale  roughness  is  due  to  the 
small  pieces  of  rubble  of  which  the  keel  is  composed. 

2.2  Intervening  ice  roughness.  Between  the  keels,  it  is  assumed 
the  intervening  ice  is  perfectly  flat.  One  consequence  of  this  is 
that  backscattering  only  occurs  at  the  ridge  keels.  This  is  consis¬ 
tent  with  observations  by  Berkson,  et_  aJL 10,  using  a  narrow-band 
scanning  sonar.  Their  results  indicate:  "(1)  very  high  level  back- 
scattering  from  well-defined  under-ice  ridges  and  (2)  very  low  levels 
of  backscattering  between  ridges."  The  roughness  of  ice  lying  between 
the  ridges  is  not  negligible,  but  it  produces  a  negligible  contribu¬ 
tion  to  acoustic  backscatter  at  small  grazing  angles.  The  authors 
speculate  that  the  coherence  of  the  forward  propagating  energy  may 
depend  on  the  roughness  of  the  intervening,  "undeformed"  ice,  since 
high  frequency  energy  striking  the  sloping  side  of  a  ridge  will 
generally  be  scattered  out  of  the  propagation  channel  at  high  angles. 


2.3  Roughness  spectrum.  All  calculations  of  scattering  strength 
in  this  article  are  based  on  the  one-dimensional  folded  roughness 
spectrum  published  by  Mellen5.  A  good  fit  to  this  spectrum  is  given 
by  the  function11: 


G1  <kx>  * 


T^x7) 


3/2 


where 

B  =  .013, 

kQ  *  .05  (inverse  meters). 

The  isotropic  two-dimensional  spectrum,  F2 ,  associated  with  is 
obtained  using  the  Hankel  transform.  This  integral  can  be  carried  out 
in  closed  form  to  obtain: 


f2  (|k|) 


_ B 

*(“!*  |k|2)2' 


where  k  is  the  two-dimensional  wave  vector: 


k  -  (kx,  ky). 


2.4  Spectrum  correction  factor.  The  roughness  spectrum  is  used 
to  characterize  the  roughness  on  the  ridge  keels.  Mellen's  spectrum 
for  the  entire  surface  can  be  thought  of  as  a  weighted  average  of  the 
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roughness  spectra  for  the  ridged  and  smooth  portions  of  the  surface. 
A  correction  factor  must  be  applied  to  Mellen's  spectrum  to  account 
for  the  fact  that  all  roughness  is  localized  in  the  keels.  This 
factor  is 

1/Pr, 

where  Pr  is  the  fraction  of  the  surface  area  that  is  ridged. 

2.5  Ridge  orientation.  The  orientation  of  the  ridges  is  assumed 
to  be  random  in  this  model.  However,  sound  propagating  in  a  given 
direction  will  tend  to  collide  most  often  with  those  ridges  lying 
perpendicular  to  its  path.  The  distribution  of  ridge  orientation  with 
respect  to  a  ray  path  is  a  variation  on  the  Buffon  needle  problem, 
discussed  in  Feller.^  Let  $  be  the  angle  between  the  ridge  axis  and 
the  perpendicular  to  the  track  of  the  ray  along  the  ice.  The  proba¬ 
bility  density  for  the  angle  4,  given  that  the  ray  collided  with  the 
ridge,  is: 

cos  #  d$,  where  jMt  <  $  <  ji  . 

2  2  2 

2.6  Shadowing  by  ridge  keels.  Shadowing  by  ridge  keels  works  to 
increase  scattering  strength  at  Tow  grazing  angles,  since  it  is  the 
flat  portion  of  the  ice  surface  which  will  be  shadowed.  The  probabil¬ 
ity  of  scattering  from  a  ridge  on  a  given  bounce  is  thus  enhanced  at 
low  grazing  angles.  The  probability  that  a  ray  will  be  scattered  from 
the  front  of  a  ridge  is  estimated  to  be: 

p  _  (h/tany  +  (n/2)  wQ)/S, 

S  with  maximum  value  1, 
where  h  is  the  average  keel  depth, 

wQ  is  half  the  average  keel  width, 

S  is  the  average  keel  spacing,  and 
y  is  the  grazing  angle  of  the  incident  ray. 

The  formula  for  P5  is  the  ratio  of  the  length  of  the  ridge  shadow  to 
the  mean  spacing  between  ridges.  The  length  of  the  shadow  is  the  mean 
distance  from  the  front  of  the  ridge  to  the  top  of  the  ridge  along  the 
track,  including  the  effects  of  orientation,  (x/2)  wQ,  plus  the 
distance  from  the  peak  of  the  ridge  back  to  the  shadow  image  of  the 
peak,  h/tany. 

2.7  Mean  backscatter ing  strength.  The  mean  backscattering 
strength  o l  a  small  patch  of  the  surface  is  the  local  backscattering 
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strength  from  a  ridge  at  a  given  orientation,  averaged  over  the  ridge 
orientation,  and  multiplied  by  the  probability  that  the  patch  is  on  a 
ridge.  Following  Bass  and  Fuks7,  the  local  backscattering  strength  of 
a  patch  on  a  ridge  at  a  given  orientation  $  is: 

S($)  ■  4  k4  (a»n)**  i—  F,  (2k  (l-(a»n)2)1/2), 

PR 

where:  a  is  the  unit  incoming  wave  vector  (cosy,  0,  siny), 

y  is  the  grazing  angle  of  the  ray, 

n  is  the  unit  normal  to  the  ridge  surface  (cos$  *  sin9, 
sin$  *  sin9,  cos9), 

9  is  the  base  angle  of  the  triangular  ridge  cross  section, 
k  is  the  acoustic  wavenumber, 

PR  is  the  percentage  of  the  surface  covered  by  ridges,  and 
F2  is  the  two-dimensional  isotropic  roughness  spectrum. 

The  mean  backscattering  strength  is: 

s  «  Po  J*/2  sC <t>>  d$. 

3  -it/ 2  2 

3.  DISCUSSION  OF  RESULTS 

3.1  Backscatter  data.  Estimates  of  acoustic  backscatter  based 
on  composite  surface  theories  are  compared  with  measured  data  from 
three  locations**  in  Figures  1  and  2.  The  data  were  taken  in  the 
frequency  band  from  1.28-2.56  kHz.  Above  320  Hz  there  is  very  little 
frequency  dependence, 

3.2  Gaussian  slope  model-data  comparison.  In  Figure  1,  the  data 
are  compared  with  the  prediction  at  1". SI  kHz  of  a  composite  surface 
model  with  Gaussian  slope  distribution.  The  prediction  lies  well 
below  the  measured  values.  Furthermore,  the  model  exhibits 
substantial  frequency  dependence  in  our  implementation. 

3.3  Triangular  ridge  model-data  comparisons.  In  Figure  2  the 
data  is  compared  with  two  predictions  at  1  .Si  kHz,  based  on  a 
triangular  ridge  model.  The  predictions,  using  triangular  ridge 
shapes  with  base  angles  of  30*  and  35*,  respectively,  give  good  agree¬ 
ment  with  observed  backscatter  levels.  Furthermore,  the  prediction  is 
independent  of  frequency  for  higher  frequencies,  in  quantitative 
agreement  with  data  in  Brown  and  Milne1*  for  five  octave  bands  from  452 
Hz  to  7.24  kHz.  The  characteristics  of  the  triangular  ridges  are: 
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Figure  1 

Estimated  backscattering  strength  (-•-•-)  at 
1.81-kHz  using  a  composite  surface  model  with 
Gaussian  slope  distribution  is  compared  with 
data  from  Brown  and  Milne4  for  Spring  and  Summer 
pack  ice  at  three  locations  in  the  frequency 
band  1.28  to  2.56  kHz  . 


OftMINO  ANGLE 


Figure  2 


Estimated  backscattering  strength  at  1.81-kHz 
using  a  composite  surface  model  with  triangular 
ridge  shapes  having  a  base  angle  of  30*  (-•-•-) 

and  35*  ( - )  are  compared  with  data  from  Brown 

and  Milne1*  for  Spring  and  Summer  pack  ice  at  three 
locations  in  the  frequency  band  1.28  to  2.56  kHz  . 
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3.3.1  Mean  height  (h)  “4.3m 

3.3.2  Mean  spacing  (S)  “  100  m 

3.3.3  Base  angle  slope  “  30*  or  35*. 

4.  CONCLUSIONS 

4.1  Diffraction  scattering  on  steep  slopes.  The  composite  sur¬ 
face  scattering  theory  based  on  triangular  ridge  keels  gives  reason¬ 
able  agreement  with  data.  The  technique  is  based  on  discrete  ridge 
statistics  and  roughness  spectra  which  can  be  measured  from  a  sub¬ 
marine  platform.  A  base  angle,  sometimes  called  the  angle  of  repose, 
of  the  triangular  ridges  of  about  35*  seems  to  give  good  results.  It 
is  the  fact  that  scattering  takes  place  primarily  from  steeply  sloped 
surfaces,  in  the  triangular  ridge  model,  which  contributes  most  to  the 
improved  agreement  over  the  Gaussian  model.  The  actual  scattering 
mechanism  is  diffractive  scattering  from  small-scale  roughness,  as 
estimated  from  the  roughness  spectrum. 
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